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Quantum effects from particle production on background 
evolution and Cardy-Verlinde formula in f(R) gravity 
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Abstract 

We investigate particle production in an expanding universe under the assumption that the La- 



o 

grangian contains the Einstein term R plus a modified gravity term of the form R a , where a is a 

o ' 
oo 



constant. Dark fluid is considered as the main content of the universe and the big rip singularity 
appears. Quantum effects due to particle creation is analysed near the singularity and we find that for 
a £]|, 1[, quantum effects are dominant and the big rip may be avoided whereas for a £ J the dark 
fluid is dominant and the singularity remains. The Cardy-Verlinde formula is also introduced and its 
, equivalence with the total entropy of the universe is checked. It is found that this can always occur 

in Einstein gravity while in f(R) gravity, it holds only for a = £^+6 ' n being the space dimension, 
corresponding to the situation in which the big rip cannot be avoided. 
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1 Introduction 

General relativity (GR) is widely accepted as a fundamental theory to describe the geometric prop- 
erties of spacetime pp. In an homogeneous and isotropic spacetime, the so-called Friedmann- Robertson- 
Walker (FRW) model, the Einstein field equations give rise to the Friedmann equations that describe 
the evolution of the universe. It describes the universe from around one second after the big bang to 
the present matter dominated era. This standard model's success is in part due to several of its predic- 
tions having been verified by observation. For example, the abundances Helium with respect to other 
light elements observed in the universe agrees well with the predictions of this model. The period of 
recombination is strongly supported by the CMB which is arguably the strongest evidence supporting 
the standard model. A good amount of observational data indicate quite clearly that the present uni- 
verse is in an accelerated expanding phase [U [31 H]. The universe may enter in a so-called phantom 
era with an effective equation of state parameter uj less than — 1. The simplest way to explain this 
phantom dark energy era is based on the introduction of a scalar field with negative kinetic energy 
El H HI El HS1 HH US US HH HSl H3 ED- The main property of such a phantom field in the 
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accelerating FRW universe is the appearance of finite-time future singularity p~9l [20] (see [2TJ [22] [23] for 
the classification of singularities). In turn, this can lead to bizarre consequences such as negative entropy 
of the universe [24] . 

However, GR is not the only relativistic theory of gravity. In the last decades several generalizations of 
Einstein field equations has been proposed [25]|26]- Within these extended theories of gravity nowadays, 
a subclass known as f(R) theory are an alternative for classical problems, as the accelerated expansion of 
the universe, instead of the dark energy and quintessence models 27, 28, [29] . Then, an effective phantom 
phase can be realized without a scalar phantom. f(R) theories of gravity are basically extensions of the 
Einstein-Hilbert action with an arbitrary function of Ricci scalar R [30] [31] [32] [33] [34] [35] . There are 
three formalisms in deriving field equations from the action in f(R) gravity. The first is the standard 
metric formalism in which the field equations are derived by the variation of the action with respect to 
the metric tensor g M „. In this formalism the affine connection r*, depends on g^. The second is the 
Palatini formalism [3 6) in which the metric and the connection are treated as independent variables when 
the action is varied. The third is the metric-affinc formalism is which the gravitational action is a general 
function of the scalar curvature while the matter action is allowed to depend also on the connection 
[51] . Note that these approaches give different field equations in f(R) gravity, while for the GR action 
they are identical with each other. The methodology leads to a boundary contribution which is usually 
dropped out setting null fluxes through Gauss-Stokes theorem [37, 38 . In this paper we will use the 
metric formalism for obtaining the field equations. 

As we mentioned above, in Einstein gravity, when the dark energy is introduced for explaining the late 
acceleration of the universe, finite time singularities can appear. Note that even in f(R) gravity where we 
do not need to introduce the dark energy for explaining the late acceleration, finite time singularities can 
also appear in a background driven by a dark fluid. The question to be asked is: does particle creation 
can avoid these singularities or modify their nature? Particle production phenomenon in f(R) is then 
analysed in this paper and its impact as quantum effects is checked at singularity time. 

Particle production from the vacuum by the gravitational field using quantum field theory in an 
expanding universe has been introduced firstly by Parker in the 1960s [39]. One of the interesting results of 
this work is that in a radiative dominated expanding universe there is no massless particle creation due to 
the conformal invariance of the metric. Latter, quantum process of particle production has been studied by 
several authors, done in the course of the cosmological expansion [40] [4T1 [42| |43 | [44 | l45l l46| [47l l48l l49| 150] . 
Recently, various investigations in the aim of dealing with singularities have been done. Batista and 
collaborators [51] studied the effects of particle creation when a massless scalar field is minimally coupled 
with the Einstein gravity. They found that the energy density of created particles never dominates the 
phantom energy density. In the same way, quantum effects near the big rip are studied in |52] where they 
used the n-wave regularization for calculating the energy density of particle creation and found that, in 
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this case, it tends to infinity when the big rip is approached and becomes the dominant component of 
the universe. This means that big rip can be avoided. Pavlov |53j computed both the number density 
of created particles and the stress-energy tensor for a conformally coupled massive field for the case in 
which uj = —5/3. He found that quantum effects are not important if the field mass is much smaller than 
the Planck mass and the time left to the big rip is greater than the Planck time. Bates and Anderson [54] 
used a background field approach in which the energy densities of the quantized fields are computed in the 
background spacetime which contains the big rip singularity. They found that for fields in realistic states 
for which the energy density of the quantized fields is small compared to that of the phantom energy 
density at early times, and for spacetimes with realistic values of w, there is no evidence that quantum 
effects become large enough to significantly affect the expansion of the spacetime until the spacetime 
curvature is of the order of the Planck scale or larger, at which point the semi classical approximation 
breaks down. 

The calculations of particle production is usually done comparing the particle number at asymptot- 
ically early times, or with respect to vacuum state defined in two different frames. However, treating 
of quantizaton in curved space the main problem concerns the field theory interpretation in terms of 
particles. Note that in curved spacetime there is no Poincare group symmetry and the notion of vacuum 
becomes ambiguous. The problem may be solved using the diagonalization method of instantaneous 
Hamiltonian by a Bogoliubov transformation which leads to a finite results for the created particles num- 
ber [49] . Particles are created because the modes with positive and negative frequency of the field become 
mixed during the universe expansion. Particle production is directly connected with the curvature of the 
universe and then when the field equation is put in the form of harmonic oscillator equation without 
friction part, its effect appears in the effective mass. Hence, in the radiative dominated universe with 
massless scalar field, either of zero or non-conformal coupling, the effective mass vanishes and the modes 
are not mixed in the course of the expansion due to conformal invariant of the metric and as consequence 
there is no creation of massless particles. Note also that in the static universe where the scale factor is 
constant and the energy density of ordinary matter is null, even the field is massive and minimally or 
conformally coupled with Einstein gravity, the effective mass vanishes and once again there is no particle 
production. However, even the density of ordinary matter is null and the modified gravity is taken into 
account, this could not hold. Under some assumptions the scale in such a condition could not be constant 
and as consequence, particle could be created. 

In this paper we use the ansatz that the function f(R) = R + XR a and analyse the phenomenon of 
particle creation. In a first step we verify a known result that the particle production is the same as in 
Einstein gravity and this is seen through the scale factor behaviour. We consider the main content of the 
universe as the dark fluid and analyse particle production phenomenon. Note that with this dark fluid 
driven the background, the universe may present the big rip. Quantum effects from particle creation 
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is then checked near this singularity, comparing the classical energy density of the dark fluid with the 
renormalised energy density from the stress tensor corresponding to a massless minimally couple scalar 
field. We find that depending on the values of the parameter n, the big rip can be avoided. Then, we 



In the second step an attention is devoted to the Cardy-Verlinde formula based on inhomogeneous 
equation of state. Verlinde |55j made an interesting proposal that Cardy formula [56] in two-dimensional 
conformal field theory can be generalized to arbitrary spacetime dimensions. Verlinde further proposed 
that a closed universe has subextensive (Casimir) contribution to its energy and entropy with the Casimir 
energy conjectured to be bounded from above by the Bekenstein-Hawking energy and as consequence, 
one obtains a very deep relation between gravity and thermodynamics 57 . Within the context of the 
radiation dominated universe, such bound on the Casimir energy is shown to lead to the Hubble and 
the Bekenstein entropy bounds respectively for the strongly and the weakly self-gravitating universes. 
The generalized entropy formula, called the CV formula, is further shown to coincide with the total 
entropy of the universe coming from the Friedmann equations. These results were later generalized 
[58l 159} l60l EU [62l IE2 EH EI] . Our goal here is to analyse the equivalence between the CV formula and 
the total entropy coming from the generalized Friedmann equations in f(R) theory. Hence, we find that 
for the inhomogeneous EoS, the CV formula can always be recovered from the generalized entropy of the 
universe in Einstein gravity but only for a = m gravity. Moreover we check the behaviour of 

the CV formula terms near the singularity and find particularly that the Casimir energy decreases as the 
universe expands and vanishes near the singularity. 

The paper is organized as follows. In the second section a breve notion on particle production in 
expanding universe is addressed with the use of the diagonalization of instantaneous Hamiltonian by 
Bogoliubov transformation. In the third section we present the view of modified gravity about particle 
production and also quantum effects near the big rip are analysed. In the fourth section the equivalence 
between the total entropy of the universe and the CV formula is checked and an analysis of the quantum 
effects on CV formula terms near the singularity is also investigated. The conclusions and perspectives 
are presented in the final section. 

2 Particle production in expanding background 

Let <fi denotes a scalar field of mass m, and a(rj) the scale factor for spatially flat homogeneous and 
isotropic FRW spacetime with conformal line element 



find that for a € 



5, 1[ the big rip can be removed whereas for a G J it remains. 
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The Lagrangian density for this massive scalar field minimally coupled with the gravitational field in a 
conformal spacetime reads 



where t]^ is Minkowski metric. The corresponding field equation reads 
We can decompose the real scalar field into the modes as 



<j){r},x) = - 



d 3 k 



^ — iki 



a J (2tt) 3 / 2 
The mode functions now satisfy the equation 



Xk 



a 
a 



k 2 + m 2 a 2 - — I Xk = , 



(2) 



(3) 



(4) 



(5) 



where the prime denotes the derivative with respect to the conformal time rj. The modes \k satisfy to 
the Wronskian relation 



Xk(v)Xk'(v) ~ X*(v)x' k (v) = i 



(6) 



The system is quantized in a standard fashion by treating the field \ as an operator, imposing the 
equal-time commutation relations 



[x(v,x),tt(t],x')} = iS 3 (x-x') 



(7) 



where ir = dx/di] = is the canonical momentum. Then, the operators ar and at satisfy the usual 



commutation relations 

[«* 4,1 =$ 3 (k- k') , ap] = [at, at ] = . 
The vacuum state is then defined as the state for which 

a £ |0) = V k . 

Other states are built up from this by acting on it various combinations of creation operators at. 
The Hamiltonian corresponding to @ is 

H= \ fd 3 x {a 2 [{cp' f + (V(/») 2 ] + a 4 mV} • 



(8) 



(9) 



(10) 



Making use of (j4]) , the Hamiltonian can be re- written in terms of creation and annihilation operators and 
the mode functions as 



// - 2 I d'k 



Pk(n) a k a -k + p k 0?)4 a -fc + ( a k a l + a l a k 



(ii) 



where Pk(r]) and Qk{v) are defined by 

Pk = (--Xk+xX +"lxl , (12) 



a 



a 



II ^ I ,/ 7, 



Qk = [~-Xk+x')[--xt+Xl)+^kXkXt , (13) 
with ujI = k 2 + m 2 a 2 . 

In general, the field may be decomposed into many different complete set of modes and each of these 
has its own vacuum state. Suppose we label one such set of modes as x( 7 ?)- Then, in terms of these 
modes the field is 

\ f (fik - r +1 
<f>(v, x) = - J -^—^ e -* kx [Xk(r,)a % + xUv)^] , (14) 

and the vacuum state is defined by |0) =0 for all k. The Hamiltonian in this case is written as 

h{ri)ajfl_ % + P;{r,)a\a}_ % + Q k {rj) (a % a\ + ajag)] , (15) 

where i\ and Q k are similar expressions as in ([12")) and (TT51) respectively, replacing \k by Xk- Because 
of completeness, the two sets of modes Xk and Xk are related by the Bogolubov transformations that 
diagonalize the Hamiltonian i.e., P k {T)) — 0, satisfying the relation 

Xkiv) = lkXk(v) + PkXkiv) , (16) 

with the normalization condition |7fe(?/)| 2 — |/?fe(^)| 2 = 1, where ^y k and /3k are constants and called 
Bogolubov coefficients. One can compare the two vacuum states by noting that the number operator for 
the barred states is N = J d 3 kata^. Taking its expectation value with respect to the unbarred vacuum, 
one finds 

<0] JV- 10> = Jd 3 k\(3 k \ 2 . (17) 

Thus, the number of barred particles in the unbarred vacuum in the mode k is \/3 k \ 2 . Similarly, the 
number of unbarred particles in the barred vacuum in the mode k is |/3fc| 2 . Since the diagonalization 
imposes Pk{v) — 0, using (|12p with Xk replaced by Xk , we get 

a! _ _ , _ 

Xk +Xk = -i^kXk ■ (18) 

a 

Substituting this result in the expression of Q k (i]), we S e ^ 

Q k (r)) = 2L> k (r))\xk\ 2 ■ (19) 

Using (|11[) and (HHJ) we obtain 

(01 24(77)^(77) + 1 10> = , (20 ) 
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and consequently 



2 Qkiv) 2 

Using now (1T3|) and (fT9|) one can rewrite (j2~Tj) as 

2 



l/^)| 2 = ~M4 ■ (2D 



lftfo)| a = 7 



(i) + " 2 M \xk\ 2 - a -ix k xi' + Axl) + W k ? i 

(22) 



4 ^Wl»| 2 2 ' 

There are two important point to be understood with respect to vacuum states in curved space. The 
first is that, in general, it is uncertain what criteria should be used in choosing the vacuum state. The 
problem is that many of criteria used for Minkowski space as Lorentz invariance and positive frequency 
with respect to a timelike Killing vector no longer apply in curved space. The second one is that, if 
there is some "natural" choice of vacuum when the spacetime begins, it does not in general correspond 
to the natural choice of vacuum when the spacetime ends. That is, the "in" vacuum state and the "out" 
vacuum state are different. This leads to particle production as can be seen from (1221) . 

Now, to calculate \Pk(v)\ 2 we have to find Xkijf) from ([5]) and Xkiv) from (fT8|). but this will be done 
explicitly in the next section where we will analyse the particle creation phenomenon when the scalar 
field is minimally coupled with the modified gravity. 

With the metric (QJ the Einstein equation leads to Friedmann ones as 

—a- = -^pa , (23) 



3 

a" (a') 2 2 2 

2 — = — k pa (24) 



a 

,2 



where p and p are the energy density and the pressure respectively, and k — 8nG. We suppose that they 
are related by a barotropic equation of state such that 

p = ujp , (25) 

with uj the barotropic parameter. Combining Eqs. (|2"3"|) . (|2"4")l and (|2"5j) . one obtains for the scale factor 
a(rj) oc r) 1 + 3oj for the conformal time and a(t) oc t 3l - 1 +^> for the cosmic time. 

Our goal is to analyse particle creation phenomenon in f(R) gravity and check its impact as quantum 
effects on a possible appearance of singularities in the classical background when this is essentially driven 
by a dark fluid. 

3 Particle creation aspect in f(R) gravity 

In this case the curvature R in Einstein Hilbert action is replaced by a function f(R) and the total 
action reads 

S =^J d A x^d \f(R) + C( 9flV) cfy] , (26) 
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where G is the gravitational constant, g the determinant of the metric <? M „ and C the matter Lagrangian 
which depends on the metric and the matter field <f>. The field equation can be derived varying the action 
(f26|) with respect to g^ v , the so called metric formalism [66] . 
The variation of the determinant is always: 



and for the Ricci scalar, one has 



SR 



(27) 



R^Sg 



(28) 



Now, since ST^ is actually the difference of two connections, it should transform as a tensor. Therefore, 



one can write it as 

and substituting in the equation (|28[) one finds 

SR = R pv 8g» v + g^USgT - W^Jg^ . 

Hence, the variation in the action reads: 
1 



(29) 



(30) 



SS = 



2k 2 

1 
2^2 

1 

2^2 



[5f{R)^—g + f{R)5^]d i x 



F{R)8R^—g- -V^gg^Sg^f(R) 



d A x 



F{R)[R^SgT + g^g^ ~ V M V^<?H - -g^SgT f{R) 



d 4 x (31) 



where F(R) = . Doing integration by parts on the second and third terms of (|3Tj) . one gets 



SS 



J V^gSiT {F(R)Rp, - \g^f{R) + [g^U - V M V„] F{R)\ d 4 



(32) 



Demanding that action remains invariant under variation of the metric, i.e SS = 0, one obtains the field 
equations as 



F(R)R^ - ^f(R)g^ - y^,F(R)+g^DF(R) = n 2 T^ V 



(33) 



where T^ v is the energy momentum tensor of the matter fields defined by the variational derivative of C 
with respect to the metric g^ v : 

2 8{y/=gC) 



-g 



(34) 



This satisfies the continuity equation 

= , (35) 

The trace of (|33|) gives 

3DF(R) + F(R)R - 2f(R) = k 2 T , (36) 

where T = g^T^ and DF = [l/y/=g) (yf=gg» v d v F). 

Einstein gravity, without the cosmology constant, corresponds to f(R) = R and F(R) = 1, so that 
the term DF(R) in ([3"6")l vanishes. In this case one has R = —k 2 T and hence the Ricci scalar is directly 
determined by the matter. In modified gravity the term \3F(R) does not vanish in (|36p . The time-time 
and space-space components of the field equation (|33[) give respectively 



3(w + l) 



3FH 2 - - (FR - f)+ 3HF = K 2 p , (37) 
-2HF — F — 2FH - 3FH 2 + i (FR - /) = k 2 P . (38) 



HF - 2FH - F . (39) 



Making use of (|25|) and combining ([37]) and (|38|) one obtains 

1 



Fi/ 2 + HF--(FR- f) 



The difference with the modified gravity is that the usual matter energy density is not the effective one. 
Then, even when there is no energy density for usual matter, there exists a kind of fluids called dark 
fluids [55] . This appears explicitly when we put ([57)) and ([3"5|) in the followings forms, neglecting the 
contributions of any other kind of usual matter, 

3i/2 = | (\( FR -f)- 3Hp ) ( 4 °) 
-3H 2 -2H=±(F + 2HF-±(FR-f^j (41) 

Comparing Eqs. (|40l) and (|4"Tj) with the standar Friedmann equations (3-ff 2 = K 2 p and —3H 2 — 2H = K 2 p), 
the both right sides of these equations may be identified with the energy and pressure of a perfect fluid, 
in such a way that the barotropic parameter of the EoS for this dark fluid is defined after simplification 

by 

. _Pdf _ -l( F R- f) + F + 2HF 



Pdf \(FR-f)-3HF 
and the corresponding EoS may be written as follows: 



(42) 



Pdf = -Pdf + ^(-HF + P) (43) 

where the subscript " df" means " dark fluid" . This inhomogeneous EoS (|4"3"|) for this dark fluid takes the 
form of the kind of dark fluids studying in several works [701 [7T1 [721 1731 174] 
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One can put the function f(R) in the Einstein- Hilbert part and a function R a as f(R) = R + XR a , 
which includes the Starobinsky's model [68] as a specific case (a = 2), A being a constant. Let us choose 
the scale factor as a(t) — a^V , where r may be determined in terms of the barotropic parameter u) as 
usually done in the Einstein gravity. Then, using the curvature as R = 12H 2 + QH, (|39j) takes a new 
form as 



1 



12r 2 + 42r J - 36r 4 - 18r\2r - l)w 



6r(2r- 1)< 2 

/I9r 2 -fir\ a „ 

. (44) 



T2r 2 - 6r\ 

+A ( - ) t 2 (-2a + 6a 2 - 4a 3 + 3r - bar + 4a 2 r - 6r 2 + 3ar 2 ) 

12r 2 — 6r ' 

-Xuj ( ^ I t z (3r — 9ar + Qa z r — Qr z + 3ar z 



t 2 

Note that in the special case of a = 1 and setting a = 1, one obtains the result of Einstein gravity, 
that is r — w^f^\ ■ Even for any a > 1, this results remains the same with the ansatz / = R + XR a . 
Remember that particle production is directly connected with the curvature of the universe. Then, since 
the scale factor remains the same as is Einstein gravity, it is obvious that the curvature of the universe 
may remain the same and also as the effective mass. This means that particle creation aspect is the same 
as in Einstein gravity when ordinary matter is taken into account. This result does not agree with one 
found by Pereira and collaborators [69]. They used Palatini formalism in the context of f(R) gravity for 
analysing cosmological particle creation for a spatially flat universe and found that a conformal invariant 
metric does not forbid the creation of massless during the radiative era of the universe. Their result is 
that, in the context of modified gravity, the scale factor in radiative universe behaves as whose of de-Sitter 
universe in Einstein gravity which would lead to particle creation and consequence implies that Parker's 
result is valid only in the context of general relativity. Note that in a radiative universe, the trace of 
the energy momentum vanishes and so, their equation (6) must lead to f(R) oc R 2 which means that 
their equation (7) would never hold. From our result in this paper, it is clear that the scale factor of a 
radiative universe in the context of f(R) gravity remains the same as in Einstein gravity. Then, without 
any approximation, Parker's result must also hold in f(R) gravity. 

However, it is important to mention that models of the form f(R) = R+XR a , that admit the existence 
of a viable matter dominated epoch prior to a late time acceleration requires that the parameter a belongs 
to]0,l[, [67]. 

Let us consider that the main content of the universe is a dark fluid, so, neglecting any contribution 
of ordinary matter. In the absence of ordinary matter fluid (p = 0), Eq ([37]) gives 

3 (1 + aAi?"- 1 ) H 2 = -{a - l)XR a - 3a(a - l)XHR a - 2 R . (45) 

The cosmic acceleration can be realized in the regime F = 1 + aAi? Q_1 >> 1 [BB]. Under the 
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approximation F aXR a 1 , by dividing Eq (|4T]) by 3a\R a 1 , one obtains 

H 2 



(46) 



Let us assume the power law expansion solutions for the scale factor, a(t) = aot r , under the restriction 
that r ^ 2. The particular case, r = 2, must be treated separately and will be presented later. With this 
power law expansion for the scale factor the Hubble parameter behaves as H = r/t and using (|46[) . we 
find 

_ (a- l)(2a - 1) 



2 - a 

The acceleration requires that we have r(r — 1) > 0, which means that 

1 -\/3 r .1 r A + V3 „, ,„ r 
a€ /=]-(»,— ^-[U]-,l[U]—^-,2[U]2,oo[ 

In terms of conformal time the scale factor behaves as 



a(7?) oc 77 9 , q = 



l-r 



(47) 



(48) 



(49) 



As we need an expanding universe, one can distinguish two different interval for the conformal time. 
Hence, for < r < 1 one has < rj < oo whereas for r > 1 or r < 0, one has — oo < r/ < 0_. We can 
now examine particle creation phenomenon in spatially flat spacetime for massless scalar field minimally 
coupled with the gravitational field, £ = 0. Then, the equation (01 becomes 



Xk 



-2 g(g— 1) 



(50) 



This is essentially the same equation that governs the evolution of gravitational waves in an expanding 
universe [75] . This equation admits solutions in terms of hankel functions: 



Xk(v) 



A k HW(k\r,\)+B k Hj»(k\ri\) 



v=\l/2-q\ , 



where A)~ and Bk are constants to be determined. Making use of the Wronskian relation 

47 

zH( 2 \z)d z HP(z) - zHi l \z)d z H^{z) = - , 
one obtains, imposing orthonormalization of the modes, 

\A k \ 2 ~\B k \ 2 = l . 



(51) 



(52) 



(53) 



For fixing the initial vacuum state we use the Bunch-Davies sate [5TJ [77] by the choice B k = and 
Ak = 1, then the solution (fBTj) becomes 



Xkiv) 



[ H^{k\v\) 



(54) 



fl 



Note that the solution (|54f reduces to the Minkowski one in which it does not occur particle creation 
phenomenon only for the time interval — oo < r\ < 0_. 

Now we can proceed to the calculation of |/3fc(?7)| 2 through the expression (f22j) . Since we are dealing 
with a minimally coupled massless scalar field, the expression (|22[) becomes 

i ($ + k 2 ) | Xfe | 2 - % (xkxV + x' k xl) + \x' k \ 2 1 



\Mv)\= A ,,, , 

4 k 2 \xk\ z 
Now, substituting ([53]) into (1551 and making the change T] 



(55) 



\Mv)\ 2 = -^v- 2 " 



r], one we get 



V 



k 

—2 {q 
k 2 rj 



H^ik^H^krj) + ffM^^-'iN - H^ik^H^kr,) - flJ a )(rfiN 



r(2) 



(2) 



(56) 



Note that at the early time, 7/ — ► oo, H^\krj) 



~) and using (1561) . we obtain 



|/3fc(oo)| 2 = 0. We see clearly through ([56]) that particle production can be known at any time and also 
the initial vacuum condition is correctly reproduced at early time where there is no particle production. 
Then, as the conformal time grows, particle creation becomes important. Note that as the big rip time 
(r/ — > 0_) is approached the rate of particle production diverges. The evolution of the particle production 
rate for aQg]|,l[is presented in Fig 1. We will see later that this behaviour of particle production rate 
has an effective impact on the singularity. 

An important problem to be put out here is that, as the conformal time 77 approaches 0_, the scale 
factor 0(77), the energy density pdf and the pressure pdf of the dark fluid, diverge: this corresponds to the 
big rip singularity. 

Obviously, the question to be asked when there is appearance of singularity is : does this particle 
production, as quantum effects, can lead to the avoidance of such a singularity? To answer to this question 
it is important to evaluate quantum energy density and pressure, and compare them with the classical 
ones of the dark fluid. The suitable expression for the quantum energy density is the renormalized one 
and this in perfectly known in the literature. For a minimally coupled scalar field, the renormalized 
energy density reads [77] 



( R 



1152tt 2 V 6 ™ 2 
1 



11527T 2 

, . - 168J2.^ + 2880Rg^ + 24R^R° 



69120tt 



-l2R a PR aP g^ - 6ARR^ + 63R 2 g^ 



Rg, 



192tt 2 C77 2 



(57) 
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Since we are leading with a massless scalar field, we need to take a null mass in ([57]) . Note that in this 
massless limit, we have to introduce an arbitrary length (or inverse mass) p 

Wff^ In ( JL) = In (Rp 2 ) - ^H^ In (6mV) . (58) 

In a full gravitational dynamical theory, there would be a term ^H^ v on the left-hand side of the 
gravitational held equations which arises from the presence on an R 2 term in the generalised gravitational 
action. Then, the final term of (1551) is proportional to ^H^ v and so may be taken over the left-hand 
side of the held equations and absorbed in the renormalisation of the coupling constant of this term. The 
other term on the right side remains in (T /il/ ) ren and is essential to the conservation of that quantity. 
Then, the energy momentum renormalised tensor of a massless minimally coupled scalar field is 



(i)h- (1)h 



11527T 2 v ^ ' 11527T 2 

1 



i> ( I + u ) +^ 



69120tt 2 



2 J \2 
- 168i? ;M „ + 288DRg„ v + 24R M(T i£ 



-12R a ^R a09flv - URR^ + 63R 2 g^j - j^g^ ■ ( 59 ) 

With the scale factor (|45p we see from (j59")l that in FRW the renormalised energy density behaves as (the 
same also occurs for the pressure) 

p ren ex rr 4(9+1) (60) 
For other hand, the energy density of the dark fluid is 



Pdf 



K 2 F 

The ratio of the energy densities is 



' ' ' ■ ocry" 2 "^ 1 ) . (61) 



^ (FR — /) — 3HF 



Pl^L r/ 2( Q -2) (g +l) (62) 
Pdf 

Since we are leading with an expanding universe we have a £ I and two important cases can be distin- 
guished: 

• For a e J =} - oo, ^^[U } 1± ^, 2[ U ]2, oo[, the product (a - 2)(q + 1) is positive and when the 
big rip time is approached , that is rj — > 0_, the ratio of energy densities goes to zero. This means that 
the dark fluid is dominant and then the big rip cannot be avoided. 

• For a £ ] h, 1[, the product (a — 2)(q + 1) is negative and the ratio of energy densities diverges. This 
indicates that quantum effects may be dominant as the singularity is approached and the big rip may 
be avoided. Note that this result agrees with that obtained in the [52] where it has been considered an 
expanding universe dominated by a dark energy fluid in the Einstein gravity context. We also mention 



13 



here that this is a case for which there is existence of viable model f(R) = R + XR a , describing an early 
matter dominated and the late time acceleration universe. 

However, it is important to mention that, the real aspect of the avoidance can be observed far before 
the singularity time. So, the question to be asked is: how long time before the big rip the quantum 
mechanism becomes important enough to begin dominating the matter content of the universe? It is a 
question that requires an analysis of many details of the evolution of the universe. But, we can answer 
this question in a semi-quantitative way. Let us put the general energy density due to quantum effects 
in the form 

Pren = PrenOV' 4 ^^ ) (63) 

where p ren o is the energy density due to quantum effects today. On the other hand, the energy density 
of the dark fluid is given by pdf = Pd/o 7 7~ 2 "^ +1 ' ) - So, their ratio is 

/ « \ 2 ( Q -2)(9+l) 

Pren PrenO I i\_ j (64) 

pdf Pdfo \r)o) 

where rjo characterizes the present conformal time. The scale factor is normalized so that it equals 1 at 
present. Let us suppose that the present-time ratio between the two densities cannot exceed the ratio of 
the total radiation observed today with respect to the total density, i.e., 10 -5 . Hence the energy densities 
due to quantum effects and the dark fluid become comparable when 

— - 10^-4(3+1) . (65) 
Vo 

This may allow estimation of how much the scale factor has increased from today till the equality moment. 
Normalizing the present scale factor to unity and re-expressing it in terms of the cosmic time, we have 

(c-l)(2 Q -l) 

where to is the present cosmic time and t s the singularity time. Using the estimates made before, we find 
that the energy densities equal at 

t s -t e ~10-^)(t s -t ) , (67) 

where t e is the equality time. Hence, typically, the energy density due to quantum effects begins to 
dominate a fraction of the present age of the universe before the singularity. We conclude that it is 
sufficient to lead to big rip avoidance. So, quantum effects prevent the divergence which would occur 
about the energy density and pressure of the dark fluid as the singularity is approached and make them 
finite. This can be interpreted as a change of dynamical regime. 

As we told above, the case a = 2, the Starobinsky's model, must be analysed separately. Then, using 
directly the equation (|46[) without any ansatz about the explicit expression of the scale factor, one obtains 

2HH + 6HH 2 - H 2 = , (68) 
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whose general solution cannot be found analytically. However, one can observe two particular solutions, 
H = and H = const, which can lead to a constant scale factor. With a constant Hubble parameter, 
the energy density and the pressure of the dark fluid are finite and then, there is no finite time future 
singularity. This result agrees with that obtained by Bamba and collaborators where they shown that 
the model R + \R 2 naturally removes finite time singularity |76j . 

4 Cardy-Verlinde formula in f(R) context 

This section is devoted to the application of the CV formula to the dark fluid characterising the main 
contain of the universe. Note that the general approach to CV formula has already been done extensively 
for an ideal fluid in [78]. They derived CV formula- like which relates the entropy of the closed FRW 
universe to its energy, and Casimir energy, for a multicomponent coupled fluid where the generalized fluid 
obeys an inhomogeneous equation of state, both for Einstein and modified gravities. It has been explicitly 
shown as special result that when the equation of state is a linear one and the barotropic parameter uj = ^ 
(the radiative universe), the CV formula is recovered from the total entropy of the universe. Also, the 
recovery of the CV formula from Fricdmann equations has been investigated in our early work |79j with 
Einstein gravity context, where the inhomogeneous equation of state contains a variable cosmological 
terms. The interesting case to be investigated here is the recovery of CV formula from the total entropy 
of the universe with an ideal fluid with inhomogeneous equation of state. As we know that the classical 
background goes toward a big rip, an analysis will be done around the singularity about CV formula 
terms. 

Looking at the equation (1431) . it clear that it is an inhomogeneous one. Then, as has been done in 
[78], one can put it in general on the form 



p df = uj(a)p df +j{a) 



(69) 



and the conservation law for the energy, 



p df + nH(l + uj)p d f = 



(70) 



where n is the space dimension, becomes, using (|69|) 

dpdfja) n(l +u(a)) 



Pdf (a) 



= —n 



3(a) 



(71) 



da a 

The general solution of ([71]) is 



a 




(72) 



where Q is an integration constant. Identifying (|43|) with ([7T]). one gets 




(73) 
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and consequently, G(a) = 0. Making use of (|45|) and setting Q — 0, one gets 

«r<a) ^)= (/_1 ) a (2a -l) < 74 > 

Hence, the total energy in the volume V — a n behaves as E = pV oc a n+9 ^ , which is also the behaviour of 
the extensive and sub-extensive energies. Assuming the conformal invariance, one can write the extensive 
and sub-extensive (the Casimir energy) parts of the total energy as 

E E = - bl () S~^, E c = n b » S-***, (75) 

where b\ and 62, and S the total entropy of the universe. Then the total entropy of the universe is 



S : 



-^E C (2E-E, 



c) 



9(<*} + l 



(76) 



Vhfo 

Then, we see that the CV formula is recovered for n + 1 = — g(a), which means that the dimension 
of the spacetime may be |<?(a)|. Then, g(a) has to take only discrete values. Solving the equation 
g(a) + 11 + 1 = 0, one finds two value for the parameter a, 

= ol 2 = - — — . (77) 

2n + 6 

Then, it appears that CV formula may always be recovered from the total entropy of the universe in 
Einstein gravity. However, this occurs in f(R) gravity only for a = ai. On the other hand it is important 
to note that for either Einstein or f(R) gravity, if n + 1 7^ —g{a), CV formula cannot be recovered. 

Let us now analyse what happens about CV formula terms around the big rip singularity. Note that 
the CV formula is reproduced for g(a) + n = —1, then, the Casimir energy is inversely proportional to 
the scale factor. We find then an interesting result that, in an expanding universe, the Casimir energy 
decreases and as the singularity (big rip) time is approached, the scale factor diverges and consequently 
the Casimir energy vanishes. This result is in agreement with the Brevik and collaborators's result where 
they shown that either with viscous or non-viscous fluid the Casimir effect fades away near the big rip 
|80) . On the other hand, when we put ot-2 on the form 2 (n+i)+A > ^ appears evidently that it is less that 
1/2 and then must belong to J. Hence the CV formula is reproduced only when the big rip cannot be 
avoided. 



5 Conclusion 

We studied particle creation aspect in f(R) gravity in which we considered the assumption f(R) = 
R + XR n . Note that particle creation is directly connected with the curvature of the universe. When the 
curvature vanishes, the modes are not affected by gravitational field and particle creation phenomenon 
cannot hold. This aspect is realized when we considered a massless scalar field in a radiative universe. 
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Also in a static universe where the scalar remains constant at any time, the curvature vanished and 
particle creation cannot be expected. Note that these situation appears clearly when the field equation 
is put on the form of an harmonic oscillator equation. The frequency depends on the wave number k 
and on the effective mass. Then, the gravitational field effect through the curvature is incorporated in 
the effective mass. Consequently, when the curvature grows, the effective mass does not vanish and the 
particle production can hold, at least when the mass is considered and/or the coupling is not conformal. 
It appears that particle creation consequently is linked with the behaviour of the scale factor. We found 
then that the scale factor does not change in the f(R) gravity with respect to his form in Einstein gravity 
and consequently the particle production aspect is the same as in Einstein gravity. We considered that 
the main content of the universe is a dark fluid and then the universe may evolve toward a finite time 
singularity, the big rip. We then analysed quantum effects near this singularity comparing the classical 
energy density and the pressure of the dark fluid with the renormalized energy density and the pressure 
due to quantum effects. Hence, we find that for n 1[ the ratio of the renormalized energy density due 
to quantum effects to the classical energy density of the dark energy diverges as the big rip is approached. 
We conclude that quantum effect are dominant near the singularity and the big rip may be avoided. For 
n£ J, the ratio goes do zero as the big rip is approached. This means that the dark fluid is the dominant 
component of the universe near the singularity and we conclude in this case that the big rip cannot be 
avoided. The Starobinsky's model has been analysed separately and reveals the absence of singularity. 

However, due to the fact that power law expansion for the scale factor presents restriction on Starobin- 
sky model, it would be interesting to use an exponential solution and incorporating conformal anomaly 
as quantum effects around the big rip. Also, in the case for which the singularity is not avoided, it would 
be interesting to allow the cosmic fluid to possess viscosity for analysing the possible avoidance of this 
singularity. We propose to address these investigations in a future work. 

Another interesting point to which we devoted our attention in this work is the equivalence between the 
CV formula and the generalized total entropy of the universe, coming from the inhomogeneous equation 
of state that appears in the f(R) gravity. We found that the CV formula can always be recovered in 
Einstein gravity while in f(R) gravity this equivalence holds only for a = £^+5 corres P on ding to the 
situation in which the big rip cannot be avoided. We also analysed the behaviour of CV formula terms 
and found that as the universe expands, the Casimir energy decreases and fades away as the big rip is 
approached. 
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Figure ll The evolution of particle creation rate as function of conformal time with the mode k — 1 and a — 2/3 
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